Elliptic Flow Suppression due to Hadron Mass Spectrum 
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Hadron resonance gas models provide a good description of the equation of state of quantum 
chromodynamics determined by lattice QCD calculations at temperatures T ~ 100 — 155 MeV. In 
this paper we investigate the effects of an exponentially increasing hadron mass spectrum (Hagedorn 
spectrum) on the azimuthal anisotropy of the rapidly expanding matter formed in ultrarelativistic 
heavy ion collisions. If the temperature at which the conversion from fluid degrees of freedom to 
hadrons is sufficiently close to the Hagedorn temperature, the production of Hagedorn resonances 
suppresses the differential elliptic flow of all hadron species. 

PACS numbers: 12.38.Mh, 24.10.Pa, 24.85. +p, 25.75.Dw 



Ultrarelativistic heavy ion collisions at the Relativistic 
Heavy Ion Collider (RHIC) and the Large Hadron Col- 
lider (LHC) are able to reach temperatures high enough 
to create and study the Quark-Gluon Plasma (QGP). 
Nevertheless, extracting the QGP properties from heavy 
ion collisions is still a challenge since this state of mat- 
ter is only created transiently. Experimentally, it is only 
possible to measure hadrons, leptons, and photons pro- 
duced throughout the collision with most of the hadrons 
being formed in the final stages of the collision. 

Therefore, in order to study the properties of the QGP, 
we have to model the whole heavy ion collision: from 
the formation and thermalization of the QGP to the dy- 
namics of the hadron rich medium formed at the end of 
the reaction. The theoretical modeling of ultrarelativis- 
tic heavy ion collisions has become very sophisticated in 
the last decade, achieving a level of precision without 
precedent. Currently, a state of the art description of 
the QGP formation and subsequent evolution should in- 
corporate: 1) a description of the pre-equilibrium phase 
of the QGP using gluon saturation inspired models [l|, 
2) the time evolution of the QGP using second order rel- 
ativistic dissipative fluid dynamics, 3) a conversion to a 
hadron rich phase in local thermodynamic equilibrium 
via Lattice QCD (LQCD) equations of state [2j-|6|, and 
finally 4) a description of hadronic matter using hadronic 
cascade simulations. For the current theoretical descrip- 
tion of heavy ion collisions see [7-21|. 

A simulation that includes at least items 2) and 4) is 
usually referred to as hybrid model (see, for instance, |7|), 
which couples fluid dynamics to hadronic transport simu- 
lations. These are thought to provide a more reliable de- 
scription of the hadronic matter formed at the late stages 
of the collision which would then remove uncertainties in 
the extraction of thermodynamic and transport proper- 
ties of the QGP. 

In practice, in hybrid models the transition between 
the fluid degrees of freedom to the hadronic ones is im- 
plemented via the Cooper- Frye method [22[ . Usually this 



procedure employs an isothermal space-time hypersur- 
face, with the distribution of hadrons computed in such 
a hypersurface being used as an initial condition and also 
boundary condition for the hadron cascade simulation. 
Note that in order to implement this process it is nec- 
essary to provide not only the momentum distribution 
of hadrons inside each fluid element but also the hadron 
chemistry of the system. At the temperatures when most 
hybrid models switch from fluid dynamics to transport 
theory, T sw ~ 130 — 165 MeV, the hadron chemistry is 
not fully known and is still a subject of intense investi- 
gation. 

Recently, LQCD thermodynamics [3J at temperatures 
T ~ 100 — 155 MeV has been shown to be compati- 
ble with calculations performed using Hadron Resonance 
Gas (HRG) models where the hadron density of states 
increases exponentially p(m) ~ gW 1 ) exp(m/Tf/) with 



lim 



m/Tff»l 



g(m) 



24j . The main parameter that 
characterizes the increase of the hadron density of states 
in this case is Tjj, known as the Hagedorn temperature 
[25L , which is an energy scale of the order of the QCD 
(pseudo) phase transition temperature In this pa- 

per, the hadronic states with masses larger than those 
measured in the particle data book [27j m > 2.5 GeV, 
whose existence is implied by p(m), are called Hagedorn 
resonances. These massive states were shown to affect 
chemical equilibration times (2^ - l30j . particle ratios [3l| . 
and the shear viscosity of a hadronic gas [32j, |33[ . While 
an experimental confirmation of an exponential increase 
in the number of hadron states may be challenging (3^ - 
|36| i the validity of such an exponential spectrum is moti- 
vated by the success of hadron models in describing low 
T lattice data in SU(Z) pure glue 37, 38| and also QCD 

[HS1. 

In this letter we show that the presence of Hagedorn 
resonances in the equation of state has an effect on the 
basic dynamical observables of heavy ion collisions lead- 
ing to an increase in the total hadron particle spectrum 
and to a suppression of the differential elliptic flow of all 
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hadrons species. The effect will be more significant in 
hybrid models in which the fluid degrees of freedom are 
usually converted into hadronic ones at isothermal hyper- 
surfaces with temperatures where Hagedorn resonances 
are highly populated. This introduces a new source of 
uncertainty to the determination of the viscosity effects 
in the QGP formed in ultra-relativistic heavy ion colli- 
sions. 

The thermodynamics of a resonance gas with 
Hagedorn-like resonances has been studied a long time 
ago [39l - l42j ]. The equilibrium pressure of such a system 
at temperature T is given by (we use classical statistics 
throughout this paper, for simplicity) 



P(T) 



2tt 2 



dmm p(m)K2 (m/T) , (1) 



where the integral is limited from below by a mass scale 
fn-min (taken here to be zero) and from above by M max . 
The standard HRG pressure computed using the discrete 
set of hadron states from the particle data group can be 
obtained from this continuum model by taking the ap- 
propriate discrete limit of the integral above. The effect 
of Hagedorn resonances on the thermodynamics can be 
seen in Fig. [TJ Following 0, [24| jtne trace anomaly ob- 
tained from lattice calculations 3] is compared to that 
of the HRG model and also to that found in Hagedorn 
gas models with density of states: px(m) — A\e m l TH1 
with A x = 2.84 GcV" 1 and T H1 = 0.252 GeV and 
p 2 (m) = A 2 e m / TH2 1 (to 2 + to§) 3/2 [HI where A 2 = 0.37 
GeV 2 , T H 2 = 0.167 GeV, and to = 0.5 GeV (with max- 
imum masses taken to infinity). One can see that the in- 
clusion of Hagedorn states into the resonance gas model 
allows fo r ag reement with lattice data up to T ~ 155 
MeV " 



23, 24 



43|. 



In a single component gas, an increase in the system's 
total energy leads to an increase of the individual energy 
of the particles and, consequently, to a higher temper- 
ature and larger particle number density. However, in 
the case of a gas with Hagedorn resonances, any extra 
energy given to the system is used to produce more and 
more species of heavier particles (according to the expo- 
nential spectrum) , and not to increase the energy of each 
individual particle species (or the overall system's tem- 
perature) [44j. Since the typical switching temperature 
used in realistic hybrid models is similar to the Hage- 
dorn temperature, this unusual way to redistribute en- 
ergy via the production of heavier resonances will affect 
the overall momentum anisotropy of the particles that 
are introduced in the hadronic transport codes. 

In fact, consider the standard Cooper- Frye procedure 
22| for an ideal fluid over an isothermal hypersurface S 
with temperature T sw . For simplicity, we take the perfect 
fluid approximation in this first study to better illustrate 
the effects of Hagedorn resonances. The particle distri- 
bution per unit rapidity for a given species of mass m a 




T (MeV) 

FIG. 1: Trace anomaly in the hadron resonance model. The 
black solid curve denotes the result obtained in the standard 
HRG model. The short dashed blue curve was computed 
using the model defined by the exponentially rising hadron 
mass spectrum pi(m) while the long dashed red curve was 
computed using p2{m) (both defined in the text). The data 
points correspond to Nt = 10 lattice data 



and degeneracy g a is given by 

dN a _ g a 
dypTdpTd(f> (2ir) 3 



(2) 



where pt is the particle's transverse momentum, <fi is 
the momentum azimuthal angle, u M is the fluid's 4- 
velocity, and p M is the on-shell 4-momentum of the par- 
ticle. The total particle distribution is given by sum 
over all the states produced at the freeze-out surface 
dN/dypTdpT = ^2 a dN a /dyprdpr- Approximating this 
sum by an integral over the density of states p(m), we 
then obtain that the overall particle distribution per unit 
rapidity is 



dN 



dyprdprt 



M„ 



dm p{rn) 
(2tt) 3 



The total particle spectrum is given by 



dN 



(3) 



I 



dN 



The elliptic flow coefficient per unit of 



dypTdpTd<ft ' 

rapidity of the hadrons emitted from this hypersurface 
can be computed using the standard event plane method 

M 



v 2 (pt) 



f 2 
Jo 



dN 



PT dpT dtp 



cos[2(<£- fa)] 



dN 



(4) 



where 4*2 is the event plane angle. The relevant pt 
range for hydrodynamical behavior in ultra relativis- 
tic heavy ion collisions is p™ ax < 3 GeV. Very heavy 
resonances with transverse masses itit — \/ Vr + m ' 2 
a few times larger than p™ ax contribute to make the 
overall distribution more isotropic since e _PM " A '/ Ts " J = 

e -m T rf/T 3m {l-p T ^/m T ) _ e -m 7 /T sra _ j n factj j t hag been 

known for quite some time that for heavy hadrons the 
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differential elliptic ffow generally increases slower with 
Pt in comparison to that found for light hadrons 11 , 4(| . 
Therefore, as T BW is brought closer and closer to T H , more 
of these heavy states are emitted and this "isotropiza- 
tion" mechanism induced by heavy resonances should 
lead to a suppression of the overall differential elliptic 
flow of the matter. 

Moreover, note that as we increase M max , more states 
are produced and the total number of particles increases. 
The pt spectrum is also enhanced and this effect be- 
comes more significant at high pt- This occurs because 
the exponential term e m / Tff in the density of states com- 
pensates the Boltzmann factor e - m Ti/T sw f or ver y heavy 
states and one obtains considerably more particles in the 
spectrum at high px by increasing M max (heavy particles 
should have flatter pt spectra in comparison to light par- 
ticles). Also, while the velocity field in the hydrodynamic 
calculation is not particularly sensitive to the change in 
the EOS due to Hagedorn resonances, note that conser- 
vation of energy and momentum through the isothermal 
hypersurface implies that these heavy states must be pro- 
duced when converting the fluid degrees of freedom into 
hadrons if the switching temperature is around 155 MeV. 

We tested these arguments by computing the elliptic 
flow coefficient of hadrons emitted from an isothermal 
hypersurface of temperatures, T sw — 130 and 155 MeV. 
The isothermal hypersurfaces were computed by solving 
(boost invariant) relativistic ideal fluid dynamics. We 
used a single averaged optical Glauber initial condition 
[47j to describe RHIC's 20-30% most central events at 
v/s = 200 GeV We further assumed that the ini- 

tial transverse flow of the system is zero. This will be 
sufficient to understand the effects of the Hagedorn spec- 
trum on the particle spectrum and elliptic flow, although 
event-by-event simulations would be required to investi- 
gate higher order Fourier coefficients. The equations of 
boost invariant ideal hydrodynamics are solved for this 
initial condition u sing a Smooth Particle Hydrodynamics 
(SPH) algorithm [50j. For simplicity, particle decays are 
not included and we fix the value of the energy density 
at the initialization time (1 fm/c) to match the expected 
number of direct pions for a given choice of the switching 
temperature. The equation of state used in this calcula- 
tion is the one presented in Q, whose low temperature 
behavior was shown to be compatible with a hadron reso- 
nance gas that includes a Hagedorn spectrum 23], |24|, |51 1 . 

In principle, the hadrons and resonances emitted from 
the isothermal hypersurface would rescatter and also de- 
cay, leading to further changes in the momentum distri- 
bution of hadrons and its anisotropy. This effect is not 
included in this work since our purpose is not to make 
realistic predictions for the values of flow coefficients but 
rather to estimate the effect the inclusion of Hagedorn 
resonances can have on such observables. Note, however, 
that one would expect that particle decays should not 
enhance the overall anisotropy of the expanding matter. 



In Figs. 2 and 3 we show the differential elliptic flow of 
all hadron species at isothermal hypersurfaces of temper- 
atures T sw = 0.130 GeV and T sw = 0.155 GeV, respec- 
tively. We assume that T sw = 0.155 GeV is the largest 
temperature at which one can still reliably say that the 
Hagedorn gas describes the lattice data. We used the 
previously defined density of states pi(m) and p2(m)- In 
both plots, the solid black curve corresponds to the value 
of elliptic flow in the case of an ordinary HRG without 
Hagedorn resonances while the short-dashed blue curve 
and the long-dashed red curve corresponds to the values 
of i>2 computed including Hagedorn resonances accord- 
ing to pi and pi-, respectively. The maximum mass of 
the Hagedorn resonances was taken to infinity although 
we have verified that the results already saturate when 
M max ~ 10 GeV. 

One can see in both plots that the addition of Hage- 
dorn resonances leads to a suppression of v 2 {pt) with 
respect to the HRG calculation. Also, note that the vi 
computed using pi is smaller than that computed using 
pi, which was expected since Tjj2 < Tffi- In fact, note 
that the elliptic flow suppression for T sw = 0.155 GeV is 
larger than that obtained when T sw — 0.130 GeV. This 
confirms our expectation that the production of heavy 
Hagedorn resonances leads to a suppression of elliptic 
flow if the switching temperature is sufficiently close to 
the Hagedorn temperature. We checked that the sup- 
pression does increase even further if T sw is taken to be 
0.165 GeV. Also, we verified that the spectrum becomes 
flatter due to the effect of resonances in the px range 
where the suppression of v^ipr) is more pronounced. 

The elliptic flow suppression discussed here can be of 
the same order of the typical elliptic flow reduction ob- 
tained due to the inclusion of rj/s ~ 1/Att viscous effects 
(see, for instance, [17|). It would be interesting to investi- 
gate the interplay between the suppression of elliptic flow 
induced by heavy resonance production and that coming 
from viscous hydrodynamic effects. Also, the effects of 
these heavy resonances on the dynamics of the hadronic 
phase as described by current hadronic cascade simula- 
tions remains to be studied. 

In conclusion, in this letter we showed that resonance 
production according to a Hagedorn spectrum leads to a 
significant suppression of the differential elliptic flow of 
all hadron species in ultrarelativistic heavy ion collisions 
if the switching temperature used in the conversion from 
fluid to hadronic degrees of freedom is close to the Hage- 
dorn temperature. The isotropization mechanism im- 
plied by heavy resonance production should amount to a 
reduction of the higher flow harmonics as well, which can 
be verified by extending the study done here using event- 
by-event calculations [18l |. Our results indicate that the 
inclusion of Hagedorn resonances in the description of the 
hadron rich phase formed in heavy ion collisions may be 
needed to improve the current estimates of the viscous 
effects in the QGP. 
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FIG. 2: Differential elliptic flow coefficient of all hadron 
species for 20-30% RHIC most central collisions computed 
using a single optical Glauber initial condition with freeze- 
out temperature T 3W — 0.130 GeV (no particle decays are 
included). The solid black curve corresponds to the value of 
the elliptic flow in the case of an ordinary hadron resonance 
gas model including all the states in the particle data book 
while the short-dashed blue curve and the long-dashed red 
curve corresponds to the values of v 2 computed adding Hage- 
dorn resonances that follow the density of states pi and p2, 
respectively. The Hagedorn temperatures in pi and p2 are 
T H i = 0.252 GeV and T H2 = 0.167 GeV. 
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